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Abstract 

We present a practical method for finding the canonical forms of arbitrary dimen- 
sional multipartite entangled states, either pure or mixed. By extending the technique 
developed in one of our recent works, the canonical forms for arbitrary mixed mul- 
tipartite entangled states are constructed, which may inherit local symmetries from 
the -|- 1 pure states. A systematic scheme to express the local symmetries of the 
canonical form is proposed, that provides a feasible way to verify the local unitary 
equivalence for two multipartite entangled states. 

PACS numbers: 03.67.-a, 03.65.Ud, 02.10.Xm 



1 Introduction 

Entanglement is one of the most important ingredients in quantum information science, 
gives impetus to the most extraordinary nonclassical applications, such as teleportation and 
quantum computation, etc [Ij. Now it is generally regarded that the entanglement is a key 
physical resource in realizing many quantum information tasks. Thus the quantitative and 
qualitative study of entanglement becomes more and more important. Though superficially 
* corresponding author 
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showing up with different features, usually not all these entangled states are functionally 
independent, they may be intrinsically the same while the entanglement property is con- 
cerned. Two entangled states are said to be equivalent in implementing the same quantum 
information task if they can be obtained with certainty from each other via Local Operation 
and Classical Communication (LOCC). Theoretically, this LOCC equivalent class is such 
defined that within this class any two quantum states are inter-convertible by local unitary 
(LU) operators [2j. 

The characterization of bipartite entangled states under LU can be well understood by 
using the singular value (Schmidt) decomposition. However things turn out to be much 
complicate when the multipartite case is concerned. The characterization of multipartite 
entanglement can be done by computing the local unitary invariants of the quantum states 
[3] . Two entangled states are LU equivalent if they have the same LU invariants, the relation 
between LU equivalence for n-partite pure states and the (n — l)-partite mixed states has also 
been noticed and is used in constructing the local unitary invariants |H |5] . The parameters 
in local invariants grow dramatically as the number of partite increases, and the problem 
of identifying and interpreting independent invariants becomes very complicated [6J. On 
the other hand, one can chose certain bases and put the quantum states in some canonical 
(standard) forms. Along this line, a canonical method was proposed in Ref. [7j, though it was 
only given in a set of constraints on the coefficients of the quantum state. Later, this method 
was reformulated into a compact form [8J. By introducing the standard form for multipartite 
states, Kraus proposed a general way to determine the local unitary transformation between 
two LU equivalent ra-qubit states [9J, however as the dimension increases, degeneracy emerges 
between the identical eigenvalues of the one partite reduced density matrix, the verification 
of LU equivalence becomes unpractical. 

Recently, in [1(J] we have proposed a practical method for finding the canonical form 
of pure multipartite state by using the high order singular value decompositions and local 
symmetry properties of the tensor form quantum states. In this work, we generalize this 
method to the mixed state case where the canonical form for arbitrary mixed multipartite 
state is constructed. And we also develop a systematic method to represent the local sym- 
metries in the canonical forms, this provides a feasible way to verify the LU equivalence of 
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two quantum states regardless of the degeneracy conditions. 

The structure of the paper goes as foUows. In section 2, we give a brief introduction 
to the basic technique of high order singular value decomposition which is used in the en- 
tanglement classification. In section 3, we reformulate the entanglement classification for 
multipartite pure state under LU in a more neat form, and a practical classification method 
for arbitrary multipartite mixed states is developed where the canonical form for mixed state 
is also constructed. After the complete classification of multipartite entangled states (pure 
or mixed) with their canonical forms, we develop a systematic scheme for verifying the local 
unitary symmetry between two entanglement classes in section 4. In section 5, practical 
examples of three and four qubits states are given. Finally, some concluding remarks are 
presented in section 6. 

2 LU equivalence of multipartite quantum state 

A general A^-partite entangled quantum state in dimension Ii x I2 x ■ • ■ x 1^ can be 
formulated in the following form 

h,l2,--- ,In 

I*) = Yl ^m2...i;vKl)K2)-Kjv) , (1) 

n=l,i2=lv>«jv=l 

where ^1^12.. An ^ *C ^-^^ coefficients of the quantum state in representative bases. Two 
quantum states are said to be LU equivalent if they are inter-convertible by local unitary 
operators, which can be schematically expressed as: 

N 
i 

^l'^2'"' '*iV 

= Yl ,«JV) ■ (2) 

21,12, ■■■ ,iN 

Here, the coefficients ipi^i^.^iN can also be treated as the entries of a tensor and hence the 
quantum states can be represented by high dimensional complex tensors. In the tensor form 
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of the unitary operator U^"''^ acting on the nth partite is defined as 



(3) 



For bipartite pure state, the tensor \1/ is a matrix \1/ = ["^jiij] ^ C''^^^^^ (matrices with 
complex numbers of Ii rows and I2 columns) where dimension of the Hilbert space for each 
partite is Ii and I2 separately. The singular value decomposition (SVD) of the bipartite 
state of dimensions Ji x I2 reads 



A = f/(i) . $ . [/(2) = diag{Ai, ■■■ ,Xi} , 



(4) 



where Aj > Xj > 0, < j, I = min{/i, I2}. A has the following two important properties 

1. the singular values Aj, z G {1, ■ ■ ■ , /} of matrix are uniquely defined. 

2. A is a diagonal matrix and uniquely defined (with prescribed order of the singular 
values) . 



In this case, the singular values of the quantum state \l/ readily characterize its entanglement 
properties under LU equivalence. Two bipartite quantum states are LU equivalent if and 
only if they have the same singular value decompositions. 

Here we introduce a generalization of SVD to high dimensional multipartite systems, the 
high order singular value decomposition (HOSVD), which was presented in [12]. 

Define the matrix unfolding of the tensor \E' G C^^^^ "^^ with nth index as 



(5) 



Here \E'(„) is a /„ x {In+iIn+2 ■ ■ ■ Inhh ■ ■ ■ In-i) niatrix. For example, 2x3x4 complex 
tensor unfolding with the second and third indexes have the following forms 



(2) 



(3) 



V^lll ^^211 

V'121 ^^221 

V'131 ^231 

/^lll ^^121 

V^112 ^122 

V^llS ^123 

V^114 ^^124 



^^112 
^^122 
^^132 

^131 
^132 
^133 
^^134 



^"212 
^222 
^232 

V'211 
V'212 
V'213 
V'214 



^113 
^123 
V'133 

V'221 
^"222 
^"223 
^"224 



^^213 i^llA 

i^233 ^134 

^231 \ 
^232 
^233 
i>234j 




(6) 
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For arbitrary A^-partite systems there exists a core tensor f2 for each tensor that is 

Q = [/(I) ® ^7(2) ® ■ ■ ■ ® ?7(^)^ . (7) 



Here f2 is a same order tensor in the Hilbert space Ji x /2 x ■ ■ ■ x J^v as and any 
order tensor fli„=i obtained by fixing the nth index to i, has the following property 

(fi.„=„fi.„=,)=5., (af))' , (8) 

where erf"'* is called the n-mode singular value of \E' and erf"'' > crj"'* > 0, V i < j . The 
singular value al'^^ symbolizes the Frobenius-norm erf"'' = ||ni„=j|| = ^/ {Qi^=i, Qi^=i), where 
the inner product {A,B) = En Eia ' ' " Ei^ ^n»2...iiv<i2...iiv (^^^ H^j for details of the 
technique). 

In the following we show how we get the core tensor by the LU transformation f/*^*^ , i G 
{!,••• , A^} in Eq.(I7]). Consider a quantum state \1/ with dimension of /i x /2 x ■ ■ ■ Jjv, another 
quantum state fl with the same dimension as \l/ is LU equivalent to \l/ if 

l] = t/W(g,[/(2)®...®;7W^, (9) 

where U^^\i 6 {1, ■ ■ ■ , A^} are unitary matrices. In the matrix unfolding form, Eq.Q can 
be rewritten as 

= [/W ■ Vl/(„) . (f/("+l. — l))T . (10) 

Here f/('^+i.-n-i) = f/(n+i) ^ jj{n+2) ^ . . . ^ jj^n) ^ ^ . . . ^ U^^'^\ and fi(„), ^(„) have 
the same dimensions: /„ rows and {In+i x In+2 ■ ■ ■ x /jy x Ji x ■ ■ ■ x /„-i) columns. Now 
consider the particular case where t/*-"^ is obtained from the singular value decomposition of 
matrix ^I/(n), i-e., 

■ vl/(„) . = diag{ai"\ ■ ■ ■ , ag)} = , (11) 

where [7*^"^ and V^'^"> are unitary matrix, and a^^^ > aj"'' > 0,Vz < j. Then Eq. flTUl) can be 
written as 
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It is clear that has orthogonal rows 

(a„=„a„=fc)=<5,fc [of^y . (13) 

And Eq. f|T3|) always holds if jg unitary matrix. In the similar way we can con- 

struct other local unitary matrices U^^\i E {1,2, ■■■ ,A^}. Eventually, all the f/*-*^ can be 
obtained, the core tensors of \E' can then be constructed via Eq.(l9]). 

From the construction of the core tensor, two of the important properties of HOSVD 
compared to its bipartite counterpart can be concluded: 

1. The n-mode singular values erf"'', i G {I,-- - ,In},n G {I,-- - , A^}, of \E' are uniquely 
defined. 

2. If Vn G {I,-- - , A^} the ra-mode singular values a^^^ are all distinct, then Q'^^^ = 
0{n)^(n) is also a HOSVD of \E' where 6(„) = diag{e*^i \ --- ,e^^''"}. Otherwise, let 
a["'^ > a^'^ > ■ ■ ■ > a^i^ > denote the distinct 77,-mode singular values of Q(^n) with 
respective positive multiplicities iJ,'f\ fi^2^\ ■ ■ ■ , /i^"^ where ^i^=i fJ-i^^ = In- In this case 
if is a HOSVD of ^ then 

^„) = 5('^) (14) 

is also a HOSVD of \E'. Here u^J G C^'^" ^''fe" are arbitrary unitary matrix and the 
diagonal blocks of S'*^"-' conformal to those n-mode singular values with multiplicity. 

From the second property it is clear that, unlike the bipartite case, the core tensor VL 
(HOSVD) of ^ is not uniquely defined. 

3 Classification under local unitary equivalence 

In this section we propose entanglement classification scheme by decomposing the LU 
equivalence of the quantum states into two correlated problems: the HOSVD and local 
unitary symmetries. First, we give a brief introduction to the entanglement classification 



(n) 



e 

i=l 



U 



(n) 
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of arbitrary dimensional multipartite pure states which was first proposed in \TU], then we 
generalize the method to the mixed states, by which the canonical forms for entanglement 
classes of mixed states under the LU equivalence can be constructed neatly. 



3.1 The equivalence under LU for multipartite pure states 

Due to the nonuniqueness of the core tensors, Q can not be identified as the entanglement 
classes of the quantum states. The philosophy of our scheme in [10] is that if we impose this 
nonuniqueness as a local symmetry within the core tensors themselves, then we can get the 
unique canonical forms. That is if we regard the core tensors fl and fl' which are related by 
local unitary operators as the same entanglement class then the HOSVD can be seen as the 
entanglement classification of the multipartite state \Ef. 

Suppose that the core tensor Q have kn distinct n-mode singular values cr^f^.i G {1, 2, ■ ■ ■ , A;„}, 
each with multiplicity of where jJ^i^ = In- Here we regard these multiplicities as 

the degeneracies of the singular values which corresponds to the case of nongeneric states 
of [n]. From Eq. fll4p we can infer that the local unitary symmetry which relates two core 
tensors takes the following form 

N 



s = 



(15) 



The core tensors Vt' and Vt related by this symmetry now can be written as 

n' = SVL. (16) 

Two different core tensors related by S belong to the same entanglement class, we can call 
the core tensors of \E' associated with corresponding local symmetry S the canonical form 
of ^. 

In order to see how this symmetry act on the core tensors we introduce the technique of 
vectorization of the matrix. With each matrix A = [aij] G C'^^^'^^, we can associated it with 
a vector A defined by 

A = [ail, ■ ■ ■ , o/ii, (ii2, ■ ■ ■ , o/i,2, ■ • ■ , dihi ■ ■ ■ ) '^hhV ■ i'^^) 



Two tensors \1/ and of Ji x J2 x ■ ■ ■ /^v which are related by local operators U^'^\n G 
{1, 2, ■ ■ ■ , A^}, can be expressed in the matrix unfolding form with the nth index 

= f/(") • ^(„) ■ (?7("+i) ® f/("+2) ® ■ ■ ■ ® U^^^ ® U^^^ ® ■ ■ • ® U^^-'Y , (18) 
With the convention of Eq.( fT7l) . the matrix equation Eq. (|T8l) . can be written as (see [T3] ) 

U(n+1) ^ u{n+2) ^ . . . ^ U{N) ^ ^ . . . ^ ^ f/W^^^^ = ygC ^'(„) . (19) 

This can be seen as a unitary transformation of a Ji x /2 x ■ ■ ■ x Jat vector to ^(„)- 
In case of n = we have the simple form of Eq. f|T9l) 



[/W ® ■ • • [/(^)^(^) = vec^( 



{N) ■ 



In this form, the symmetry between the core tensors can be represented as 



(20) 



TV 



(N) 



e 



u 



(n) 



(TV) 



(21) 



"=1 Un=l 

where Uj^ is unitary matrix of fi^J^ x /ij"-* and Yl'j"=i /"j""* = -^n- In the blocks diagonalized 
form, Eq.( l2T]) is 
















V 




which has been state in [10] 



(TV) 

JV 



■ ^(AT) = ^'(Af) 



(22) 



q(") 



We can set if the multiplicity ^j^^ = 1. In all, we have the following theorem 



Theorem 1 The core tensor Q associated with the local symmetry group S is the canonical 
form of the multipartite pure state and is the entanglement class under LU equivalence. 



From a more general point of view, any quantum state \E' associated with a local unitary 
group U = f^^*-* can also be seen as the unique entanglement class. In this unified scenario. 
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we can compare the symmetry U and S corresponding to \l/ and Q separately 



S 



(n) 

Jn 



N r kn 

e 

® u^'^ 



n 



32 



(N) 
7/ . 



31,32, 
N 



JiV = l 



u 



Jjin) ^ 



(23) 
(24) 



n=l 



The total matrix dimension of these two groups of 5" and U are 

N N 



Bim[S] 



31,32, ■■■ ,jN n=l n=l jn 



(n) 



AT 



J]j„ = Dim[[/] . 



(25) 



n=l 



Through this comparison we can see that, a arbitrary quantum state \1/ under LU trans- 
formation U forms an entanglement class U \1/ which is also a orbit under the transformation 
U. After the entanglement classification we arrive at the canonical form of the quantum 
state, the admissible transformation S which keeps the canonical form invariant is regarded 
as the symmetry of the canonical form. In Eq. fl23|) . S can be seen as a conjugate class of 
U under a particular local unitary transformation. In the special case that all the singular 
values are distinct for each partite, the symmetry becomes 



N 



s 



e 

n=l \Jn=l 



(26) 



which is just the conjug ate class of (g)^^^ U^'^ under (g)^ti U^'^ ■ 



From the quantum state point of view, tensor Q now is decomposed into several invariant 
subtensors (denoted it with u) of the Hilbert space of /i x /2 x ■ ■ ■ x /^v under the transfor- 
mation S. The dimensions of these subtensors conformal to the direct summed subgroups 
of S in Eqs. fl22l [23]) . For example Eq. fl22|) can be written as 



(1) 











u 



(1) 








u 



(1) 



(N) 

N ' 



( ^ri \ 



r2 



(27) 



where Qn and cJ^. are the segments of the column vectors ^{n) and with the dimension 
conformal to the diagonal blocks u^^®- ■ - ^u^^K u are just the vector forms of the subtensors 

CO. 



3.2 The equivalence under LU for multipartite mixed states 

The classification of the entanglement for mixed states is generally believed to be more 
complicated than pure states in many cases. However in the case of LU equivalence, it has 
been noticed that ri-partite pure state is related to its {n — l)-partite mixed state [5]. Here 
we generalize our entanglement classification method developed for multipartite pure states 
to the case of arbitrary dimensional A^-partite mixed states. 

Consider a mixed A^-partite quantum state p which is generally expressed as 

p = J2p'\'l^^)('l'i\ ' (28) 

i 

where ^jP^ = 1, Pi G M+, lipi) are partite pure states. We add an additional 0th partite 
to the original A^-partite mixed state p and formulate a A^ + 1 pure quantum state in the 
following form 

^o = 5Zp.N)|^.) , (29) 

i 

where \i) is the 0th partite. For this quantum state, we have the following fact 

Tro[|^o)(^o|] = I]p^P,(^K)IV^.)(^jlO» 
= ^PiPj{j\n){n\i)\iJi){iJj\ 

= 5^p?|^.)(^.|=P- (30) 

i 

And further we have, if \&o = U^'^^ ® E^^'^ (g) ■ • ■ (g) i?(^+^)\l>o = ^7^°^ "^o, where E is unit matrix. 
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then 



= p = Tro[|^o)(^o|] . (31) 
Prom the above two facts we can state that the following relation 

P = Y.P'Mr) ^ = Y.P^m^) (32) 

i=i Tro 1=1 

forms a bijection between p and ^o- Define this bijection as a map between ^-nd p, we 
have the following proposition 

Proposition 2 An arbitrary dimensional mixed N -partite state p' is LU equivalent to p, 
i.e., 

p' = [/(I) <^ <^ . . . <^ f/Wp [/(l)t <^ [/(2)t (^...^ U^N)^ (33) 

if and only if its pure state counterpart is LU equivalent to ^q, that is 

% = f/(°) [/(^) • • • ® U^^'^^o ■ (34) 

Proof: The relation between partite p and (A^ + l)-partite is 

+0 

P = I]p'l^i) ^ v]>o = ^p^|,)|^^) , (35) 

1=1 Tro i=i 

+0 

p' = Y.^M) n = E5^-i^')K)- (36) 

i=i Tro j=i 

Pirst, if 

p' = ® ® . . . ® [/(^V f/^'^^ ® f/^'^^ ® ■ ■ ■ ® u^""^^ 

i=l 

= Ep^I^O(V':I, (37) 
1=1 
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where = U^^^ ® f/^^) 



U^^^ipi), then \1/q correspond to p' is 



That is \l/o is LU equivalent to \E'o- 

Second if % = U^^'^ ■ ■ ■ ® U^^^^o, then 
p' = Tr,[\%){%\] 
= Tro 



(3^ 



U=ij=i 



Trn 



J2 PiPjm^'^ 



= [/(I) ® • ■ ■ ® [/(^V [/(i)^ ® • • ■ ® /7(^)^ . (39) 

Here, we have used the fact of Eq.f l3T|) . That is p' is LU equivalent to p. QED. 

We can now conclude that: if p' is LU equivalent to p then their corresponding pure 
states \E'q and \E'o can be related by LU operators; if is LU equivalent to then their 
reduced matrices p' and p are LU equivalent. We may construct the core tensor Qq from \E'o, 
then we trace out the 0th partite from the core tensor Qo and obtain the canonical form for 
p, that is T = Tro [|^o)(^o|] 

Theorem 3 The canonical form T is of the entanglement class of mixed state p up to a 
inherit local symmetry from Qq . 

This method provide a sample way to construct the canonical form for the mixed N- 
partite state p: first construct the partite pure state from p; then compute the 

core tensor Qq of \l'o; finally we arrive the canonical form by tracing out the 0th partite 
T = Tro|^^o)(^^o|. 
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4 The local symmetries of the canonical form 



We have constructed the canonical forms for both pure and mixed muhipartite states. In 
all, the construction of the canonical forms will result in a general form of Eq.( l22l) whether 
the state is pure or not. With this detailed form of the symmetries, in this section we develop 
a practical scheme to verify the LU equivalence of two quantum states which have the same 
singular values and same degeneracies for each partite. 



4.1 A general form of the local unitary symmetry 



(n) 

We start from a general case, that is we have A;„ distinct n-mode singular values a^J , in G 
{1,2,--- ,kn}, each with multiplicities of //j-""* where ^^"=1/^!"^ 
singular value vector a*^"^ for the matrix unfolding form of 



In- Define the rz-mode 



Mi 



(n) (n) (n) (n) 

■ ■ ■ , (Ti , (T2 ,^2 1 ' ' ' y<^2 ) 

■V ' ^ V ' 

(n) 



(n) 
M2 



{n) tT 



(40) 



(n) 



where a^^^ > af^^ > 0, V z < j. The local symmetry corresponding to this partite is 



v 



(41) 



where u[^\ in G {1, • " ' ; ^n} are unitary with the dimension of x fi'^'^\ 
The total local unitary symmetry S = ■ S"^*^ of the core tensors is 



V 



® • ■ ■ (8) 



4?/ 



V 



u 



(N) 
kjsi 



(42) 



which is just Eq. fl22p . Similarly the total singular values can be expressed as a matrix 

S^{a«,a(2)^...^^{iv)|_ (43) 
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We can call it the "singular value matrix" of the core tensor which is uniquely defined 
according to the properties of HOSVD. Quantum states with different singular value matrices 
are apparently LU inequivalent. Then the core tensors which have the same singular value 
matrix belong to the same entanglement class if and only if they satisfy Eq. fl27|) . In Eq.f l27|) . 
the verification of the LU equivalence of two core tensors turns to finding solutions for the 
following equation group with varying r 



(44) 



This can be seen as a fine grained LU classification problem of the subtensor Ur- Thus we 
can pick the sub tensor uir out of the tensor fl and do the HOSVD to it recursively. 

As the fine grained process goes, the recursion will end at two cases: 1, the singular 
values are all distinct for all the partite; 2, the singular values are all the same for all the 
partite. 

For the first case, if Vi„ G {1, • ■ ■ , kn} and Vn G {1, ■ ■ ■ , N}, the singular value multi- 
plicity fi'f^ = 1, then kn = In and 



Jjin) 



Eq.( l27|) turns to 



/exp(i0i"^) 


V 



exp[^(^« + ^;f + + 







exp{i92 




in)-. 



\ 





(45) 



(46) 



Here we write Ui-^i^.-ij^ instead of Uij^i^...ij^ because 0Ji^i^...ij^ now is a complex number in this 
case. A log operation on Eq.( l46p will make 



^lil + ^2i2 + ■ 1- ^- 



—i log 



(47) 



These are /i x /2 x ■ ■ ■ x Jjy linear equations for Ii + I2 + ■ ■ ■ + In phase variables, it can 
be verified immediately whether they have consistent solutions. The quantum states is LU 
equivalent if and only if there is at leat one solution to this linear equation group. 
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4.2 Completely degenerate state for all the partite 

In the completely degenerate state, the reduced density matrix for each partite is pro- 
portional to unit matrix. Consider a arbitrary A^-partite pure state with dimension of 
Ii X I2 X ■ ■ ■ In. The complete degenerate state is that VnG {1, 2, ■■■ , A^} 

-'n 

The core tensor has the following form 

^ = Nn) I^ni2---i„-ii„j„+i---iiv Nl^2 ■ ■ ■ "^n-l^n+l ■ ■ ■ ^Af) 

= 5^|z„)|a;(:"J),ne{l,---,iV}, (49) 

where {u^'^'"^\u^^"^) = j-di^i'^. The local symmetry 5* takes the following form 

S ■ fi(^) = (g) /7(") ■ Vl^N) = ^[n) ■ (50) 

n 

An arbitrary unitary matrix is unitarily equivalent to a diagonal matrix, that is 

where X^"^ is unitary matrix and = diagle**^! \ ' ' ' , e*"^^"'} is the conjugate class of U^^\ 
Eq. (!50!) now turns to 

I ■ (9) x(") ■ = (9) a:(") ■ n' . (52) 



The Eq. (]52il corresponds to Ii x I2 x ■■■ x homogeneous equations, which in the detailed 
form one of the of these equations of Eq. flS^ looks like 



«1---«JV 

Here we represent Xij as the elements of matrix X. This is a typical equation group of 
Ii X I2 X ■ ■ ■ X In equations for + /| + ■ ■ ■ + complex parameters x-"], (note we first 
solve the parameters xf^^ then impose the unitary condition on the matrix X^^^). 
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For this kind of nonlinear equations there exist simple tool called "linearization" or "relin- 
earization" [lll|15]. The key algorithm rely on the fact that for N > 2 multipartite quantum 
states, when the dimensional or number of partite increases, the number of equations grows 
much more quickly than the number of the parameters. Generally Eq.f l53p would turn out 
to be a over defined system of equations which mean that there are more equations than 
unknow parameters. 

The linearization technique goes as follows. Regard each monomial of the matrix elements 
as a individual variable 



(1) (2) (N) 



(54) 



then there will be (Ji x /2 x ■ ■ • x Jtv)^ such variables v. Eq.f l53p now can be written as 



. 



(55) 



For the sake of simplicity we use the convention that i represents the value of the bit 
string {iii2 ■ ■ ■ ^n), i-e., i = 1 = (11 •■■1) and i = 2 = (11 ■■■2), etc. Define Uji = 
e^'^'^'ji '''J2 '"'''jN ' — uj[ where j = {jii2 ■ ■ ■ Jn)- Eq. fl55|) can be reformulated as 



(56) 



where the dot means the summation over i. Taking a 2 x 2 x 2 system as an example, we 
have 



(57) 



There are 8 such equations for j runs from 1 to 8. The solution can be expressed as 



1 




V J 



j3 







V / 



+ ■ ■ ■ + Ci 



-j8 








V 1 / 



(5^ 
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where cji are new parameters. Clearly, Eq.( l58|) is a under defined equation group for param- 
eters z/ji. However there are additional equations between the products of ly^i s, i.e., 



This relation is inherit from Eq. (l5^ and the following fact 



(59) 



.(1) 










•Xj -f -1 




.(1) 

nji 


. 

'"nJn 


inJn 




. .^W . 


*jvJjv 


(60) 



X 

= X. 

For example in 2 x 2 x 2 system we have z^iii,iiiz^iii,i22 = '^iii, 121^^111, 112 or simply VuUn = 
J^i3i^i2- This imposes additional equation between parameters z/ji, and can also be viewed 
as an equation in the (smaller number of) parameters cji expressing them. The new system 
of equations derived from all the possible relations of the type of Eq. fl59l) . In solving the 
equations on cji we can use the linearization method recursively. 



Here we give a explicit formula for how many constrains of Eq. fl59|) there will be. As 
there are (/i x J2 x ■ ■ ■ x /^r)^ matrix elements, we can get (/i x /2 x ■ ■ ■ x /at)^ new parameters 
i/ji. If we multiply m times of Uji, that is 

yji ■ - i^yi'^ , (61) 

m 

we will have 

(62) 

different productions. On the contrary, according to the productions of x'l^\ the actual 
number of different productions is only 

N 



i=l 



which is much less than the number of equations (here = ]T(04y )- The number of Eq.(l62 
is greater than that of Eq. fl63|) when m > 1. For the case of 2 x 2 x 2 and m = 2 we have 



Cf2x2x2P+2-i = 2080 , (C2.+2-i)' = 1000 , (64) 

which means that we have 2080 different i^jiZ^j'i'S, but only 1000 are independent. A consid- 
erably large amount of constraint equations like Eq.( l59|) are obtained. 
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5 Examples of the canonical form for three and four 
qubits state 



Here we give two simple examples of how we can get the canonical forms of the arbitrary 
quantum state, and how we can verify whether two quantum states in the canonical forms 
can be related by symmetry S. As the entanglement classification of the mixed states can 
be reduced to specific pure states case, here we only give examples of pure states. 

We can randomly generate a 2 x 2 x 2 pure state ^ with its matrix unfolding 

_ /0.0260603 1.05491 -3.69051 0.43771l\ , . 

~ I 1.25266 1.07259 3.2378 1.5625 j ' ^ ^ 



From the algorithm of Eq. (1121) , the singular matrix is 



af^ af^ _ /5.03906 5.31586 5.17055\ 



0"2 ^2 

The core tensors is then 



a?^ afw V2.27534 1.5202 1.95825^ 



_ 1^-5.01792 0.2815 -0.354882 -0.0862168\ 
~ 1^ 0.19519 1.72088 -1.17941 -0.886923 J ' 

We give another example of four qubits state. Two 2x2x2x2 quantum states 

1 A000002 0' 

^(1) = 



(66) 



(67) 



^ ^0 1 2 
, _ 1 A 000002 0\ , . 

" 71 Vo 1 -2; ' ^ ^ 

which are already the core tensors. The singular value matrices for them are the same 

,(1) J2) J3) /I 4 4 1\ /^'(l) ^'(2) ^'(3) 

(69) 



0-1 Cr^ (Ti CTi \ 2 5 5 2 \ /^l ^1 ^1 



^Jl) J2) (3) (4)/ \1 1 1 1/ \ HI) /(2) /(3) /(4) 

"2 "2 "2 / \2 5 5 2/ \"2 '^2 ^2 ^2 / 

In the vector forms of the matrices unfolding of \E'(i) and ^(i), the symmetry S takes the 
following form 

/ei(?f'+*ef)f;(4)^f;{l) \ 

e^ef'+^ef/yW ® t/« 

e^'?^''+^'Vw ® [/(I) 

e^^2''+^^f ^7W®[/(i) 
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The core tensors are then divided into four segments correspondingly 



Ul = ■ 


-^diag{l, 0,0,1} 


uj'-- 


-^diag{l, 0,0,1}, 


(70) 


Ui = ■ 


-^diag{0, 0,0,0} 
V 10 


U2 = ■ 


-^diag{0, 0,0,0}, 
V 10 


(71) 


Us = ■ 


4=diag{0, 0,0,0} 
V 10 


^3 = ■ 


-^diag{0, 0,0,0}, 
V 10 


(72) 


UJ4, = ■ 


4=diag{2,0,0,2} 
V 10 


= ■ 


-i=diag{2,0,0,-2} , 


(73) 



Take all the above equations into Eg. ( 1271) we have the following effective equations 

g..f)+..f)f;(4)^^(i)^^^^. ^ (74) 
gi4')+^er^(4)^^(i)^^^^/ ^ (75) 

The components of Ui, u[ and u^, u'^ in Eqs. ( 170)1731) bring a contradiction to Eqs. ( 17^11751) . 
Now it is clear that there is no solutions for [/^^^ and U^^^ ^ and thus the four qubits states 
and are LU inequivalent. 



6 Conclusions 



In summary, by using the tensor decomposition method we generalize the entanglement 
classification under LU equivalence to arbitrary dimensional multipartite mixed states. The 
classification actually reduces to the construction of the canonical forms of the corresponding 
+ 1-partite pure states. With the analysis of the local symmetry in the canonical form, 
the core tensor can be decomposed into a series of subtensors which are transformed inde- 
pendently under the local symmetry. Base on this recognition of the entanglement structure, 
a practical scheme is also developed for the verification of local unitary equivalence of two 
multipartite entangled states. 
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